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PROPAGATION OF THE TRANSVERSE NORMAL STRESS 
IN A THICK PLATE DUE TO DISTRIBUTED 
LATERAL IMPULSIVE LOADINGS 

By Larry W. Oline* 

Langley Research Center 

SUMMARY 

The purpose of this theoretical investigation is to determine the magnitude of the 
elastic stresses produced in an infinite elastic plate when an impulsive load is applied 
normal to the plate surface. The plate is assumed to be infinite and elastic with an axi- 
symmetric normal load applied over a circular area. The first part of the analysis is 
devoted to establishing the equations of motion and the boundary conditions and trans- 
forming these equations by use of the Laplace transform. The normal stress distribu- 
tion is determined in particular. The second part of the analysis deals with the inversion 
of the transformed equations and is accomplished by use of an expansion method of 
Cagniard which eliminates contour integration difficulties. At this point in the analysis 
the various wave fronts which contribute to the normal stress can be identified and are 
recognized. 

The normal stress due to the initial, unreflected, irrotational wave is evaluated 
on the axis at the wave front, and the effect of varying the radius of the axisymmetric 
load is determined. The results compare favorably in the limit with the published 
results for normal stress due to a point load. It should be noted that only the normal 
stress at the wave front and on the axis is determined due to the mathematical complex- 
ities of the resulting inversion integrals. The solutions are obtained for both step and 
Dirac delta time variations. 

INTRODUCTION 

When a projectile traveling at relatively high velocity strikes a plate, elastic 
stress waves are usually initiated which propagate toward the back face of the plate. 

• . t 

These elastic stress waves are strong compression waves which decay in magnitude as 

*This analysis was performed under the NASA-ASEE Summer Faculty Fellowship 
Research Program. The author is now Associate Professor, Department of Structures, 
Materials, Fluids at the University of South Florida. 


they propagate. When such waves arrive at the stress-free back surface, they are 
reflected as tensile stress waves, which, depending on the intensity of the wave and its 
duration, may cause fracture to occur. This process is termed spallation and can pre- 
sent a problem in some situations. One example of a problem area is the main wall of 
a double-walled meteoroid bumper system. The outer wall breaks up the impacting 
meteor into an expanding gas cloud which then impacts the main wall over a certain dis- 
tributed area. Because of the intensity of loading and its duration, there is a possibility 
of spallation damage occurring. 

Problems related to stress-wave reflection at a free boundary have been 
encountered by seismologists, with some pioneering work done by Cagniard in refer- 
ence 1. A bibliography through 1960 of other authors having done work in this area is 
given by Miklowitz in reference 2. Several authors (Aliev in ref. 3, Kinslow in ref. 4, 
and Blake in ref. 5) have approached this problem by using a pressurized cavity as the 
source of a spherical disturbance along with an image system to preserve the stress- 
free boundary condition. Problems arise, however, in exactly matching these stress- 
free boundary conditions. 

Another approach to the stress -wave propagation and reflection problem has been 
to consider a finite thickness plate with a normal force applied to one side. Some pre- 
liminary work for a distributive force on an elastic half space was done by Huth and Cole 
(ref. 6) and by Thiruvenkatachar (ref. 7). Further analysis by Broberg (ref. 8), Davids 
(ref. 9), Pytel and Davids (ref. 10), and Rae (ref. 11) included the effects of the stress- 
free surface. All of these plate analyses considered the load to be applied at a point. 
Distributive loading on a plate was formulated by Thiruvenkatachar (ref. 12) who obtained 
the Laplace transform solution without completely finding its inverse. 

The present analysis considers a uniform normal load distributed over a circular 
area. The normal stress of the incident dilatational wave at the wave front under the 
load center is presented as a function of the radius of the loading. Numerical results 
are presented for various values of the radius of loading while maintaining a constant 
total load. A comparison is also made with the case of a point load to show the rapid 
decline in normal stress with increased loading area. 

SYMBOLS 

a radius of the load applied 


Cj dilatational wave speed, 


/ 


A + 2j it 
P 


2 


c 2 

F 


distortion wave speed, 
total force 



2 

Fjto) . (eq. (2 1)) 

(l + 2q 2 lj 2 ) - 4q 3 T) 2 (/l + l) 2 ^ + q 2 n 2 

h thickness of the plate 

. .,I 12 integrals 



P intensity of applied impulsive pressure 

p Laplace transform variable 

c 2 

q = — (nomial value of 0.5 for aluminum) 

C 1 

r,0,z cylindrical coordinates 

t time 

u r ,u z displacements in radial and normal directions, respectively 

u = cos \p 

y = -ii? 


constants 


Dirac delta 



9u r u r 9u 7 

A = — — + — + — - (Dilatation) 
9r r a z ' ’ 





Lame constants 


C 1 


p mass density 

°zz ,0 rr ,<y rz ,<y 99 stress components 

i p = Dummy variable in Poisson integral representation of Bessel function (see eq. (19)) 


fl = 




(Rotation) 


Subscript: 


P dilatation wave 


A dot over a symbol denotes differentiation with respect to time. 

ANALYSIS 

The problem considered is the determination of the stress response in a homo- 
geneous, isotropic, infinite, elastic plate subjected to an impulsive pressure. The plate, 
of thickness h, is loaded with an impulsive pressure of intensity P. The pressure is 
uniformly distributed over a circular area of radius a on one face of the plate (see 
fig. 1). A system of cylindrical coordinates (r,0,z) is defined with the origin at the 


4 





center of the circular area of loading and the positive Z-axis normal to the plate and 
directed toward its other face. Due to symmetry of the loading, the problem is inde- 
pendent of 6. The displacements in the r- and z-directions are denoted as u r 
and u z , respectively. 


Basic Equations 

The problem is governed by the following set of basic equations: 
Equations of motion (see ref. 13) 


a 2 

/. o v9A 90 0 u r 

( X + 2H) 5F +H S -P-J ? - 



where 


A = 



c^z 

9z 


( 1 ) 

( 2 ) 

(3) 


9u r 9u z 
9z 9r 


and A are Lame constants and p is mass density. 
Stress, displacement equations 


(4) 


9u r 

°rr = >a + 24 -jf (5) 

a 09 = + 2fi -p- (6) 

/9u_ \9u 7 \ 

a rz = V + -jf ) ( 7 ) 

9u_ 

a zz = AA + -fr- (®) 

Boundary conditions 

cr rz = 0 (z = 0; z = hj t = 0) (9) 

a zz = 0 (z = h; t = 0) (10) 
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a zz = PH(t)[l - H(r - aj] 

a zz = P5(t) |l - H(r - aj] 

where H(t) is the unit step function and 
Initial conditions 



(z = 0 for step function loading) (Ha) 

(z = 0 for Dirac delta loading) (lib) 

6(t) is the Dirac delta function. 


o 

All 

o' 

II 

o 

II A 
N 
II A 

cr 

(12a) 

o' 

All 

o~ 

II 

0 g z g h) 

(12b) 


Laplace Transformation of the Basic Equations 

Equations (1) to (4) were transformed in time using Laplace transforms and were 
arranged to form two partial differential equations in terms of the Laplace transform 
of A and These equations are 


and 


a 2 n l an 

ar 2 + r a r 





8r 2 r 3r C j2 az 2 


(13) 


(14) 


where p is the transform variable. These equations were solved by the method of sep- 
aration of variables and superposition (see appendix A and refs. 7 and 12). The expres- 
sions for the transform of the normal incident stress a zz can then be obtained by uti- 
lizing equations (3), (4), and (8). 

Since the type of damage during spallation is caused by a reversal of the initial 
compression wave to a tensile wave at the rear free boundary of the plate, the normal 
stress, particularly along the Z-axis, is thought to be the major contributor to this 
damage. The transformed expression for the normal stress due to a Dirac delta loading 
function, is 


zz 



M(z, £) Jj(aO J 0 (r£) d£ 

G(i) 


(15) 


The functional relationships M(z, £) and G(£) are defined in appendix A and Jj(a£), 
J 0 (r£) are Bessel functions of the first kind. The transformed solution for a step loading 
is obtained from equation (15) by replacing P with P/p. 
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The solution as given by equation (15) can be expressed as a summation of terms, 
each of which can be identified as reflected and unreflected shear and dilatation waves. 
This expansion is given as 



The values of Q n and Q n * are given in appendix B, along with the details of the 
development of this equation. As can be shown, Qq represents the initial dilatation 
wave, Qq* the initial shear wave, Q^ the dilatation wave after one reflection, and so 
forth. The result given by equation (16) has been previously obtained by Thiruvenkatachar 
in reference 12. 


Inversion of Laplace Transforms 

The present analysis considers the normal stress a zz , caused by the incident 
dilatation wave, and is represented by Qq in equation (16). Therefore, the transformed 
stress which must be inverted is 

]*“ (n^) J 0 (r » J l (a « e ‘“ Z d « < 17 > 

where Nj and N g are functional relationships found in appendix A. By introducing 

c n 

the variables 77 = and q = — and by making use of the integral representation 

P c j 

of the product of the Bessel functions (ref. 14), equation (17) becomes 



where the trigonometric variables R, x> and <P are shown in the following sketch and 
satisfy the following equations: 
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In addition, the Bessel function in equation (18) can be replaced with its Poisson integral 
representation (ref. 15) to yield 


_ aPp^ 


ZZ O o 


Re 


P (a - r cos <p) f ' sii^ f°° e ^ 

J 0 J 0 J 0 L . r>~2-9.\ 2 


- + i7jRcos \{^ 


^1 + 2q2p2^ . 4q 3 r?2^1 + rfi \/\ + q2r)2 


dr? di^ d<p (19) 


Equation (19) can be recast into a more simplified representation by introducing the var- 
iable (see ref. 11) 


c 




+ ir?R cos 4/ 


( 20 ) 


As rj varies from zero to infinity, t is represented by a curve in the complex plane 
from z/cj to °° (shown in fig. 2). It is observed that r\ = 0 defines the wave front 
z = Cjt. The path of integration can be transformed (see ref. 1) so that it lies along the 
real axis in the t-plane, and equation (19) becomes 





( 21 ) 


where 


Fjfa) = 


r?(l + 2q 2 7? 2 ) 

(l + 2q2? )2) _ + rp \J\ + q 2 7)2 


and H^t - is the unit step function. 

The well-known Laplace transform of the second derivative of a function is given 


by 


X 


{$} = p2f " " pf(0) " f,(0) 


( 22 ) 


da 

Therefore, noting that a zz (0) =0, 0 It < — , and ■ (0) = 0, the inverse of equa- 

C j dt 

tion (21) can be written as 
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R cos i \j 


Figure 2.- Variation of tj with t in the complex plane. 


= ^2 H ( t - q) Re <S X (a " r cos ** X Fl(T?) * sin ^ 


Inasmuch as equation (23) represents a stress wave propagating in the positive z-direction 
with velocity Cj, it is the incident dilatation wave. For a step loading, equation (23) 
would contain only one time derivative. 


1 


Stress at Wave Front 


If a fracture occurs during the reflection process, it is most likely to occur on the 
axis where r = 0. Setting r = 0 in equation (23) yields 



a 2 P 

7TC., 2 


Re 


air 


d?2 


dt 2J 0 F l<”> % d* 


(24) 


where the subscript P denotes dilatation wave. For this condition R = a, equation (20) 
becomes 


Cjt = z\Jl + rft + i?ja cos \f/ (25) 

It will be observed from equation (25) that as long as the value of a is nonzero, rj is 
a complex function of \p, which makes F ^(rj) and drj/dt in equation (24) also functions 
of 4/. These functional relationships must be determined before the integral can be 
evaluated. The value of a is assumed to be a positive number and a = 0 is treated 
as a special case. Solving for r} from equation (25) yields 

-iacjt cos x// ± zjc^t 2 - (z 2 + a 2 

V = — 1 

z 2 + a 2 cos 2 i// 

In this expression the positive sign must be chosen since tj = 0 defines the wave front 
z = Cjt. From equation (26) it can be seen that 77 is, in general, a complicated function 
of ;// which makes the integration shown in equation (24) a formidable task. In order to 
make the problem somewhat tractable, the solution was sought at the wave front. 

Case of a Point Load Where a = 0 

For the special case of a = 0, T) is real and independent of rp and can be taken 
outside the integral in equation (24). Equation (24) can then be written as 

( CTzz )p~ 27TCJ 2 dt 2 

where 7ra 2 P is replaced by F in the limit. 

Step load at a point . - The normal stress at the wave front 
is obtained from equation (27) by replacing the operator d 2 /dt 2 


z = Cjt for a step load 
by d/dt. The resulting 


-T7T?(l - 2q 2 7? 2 ) 

(l - 2q 2 7? 2 ) 2 + 4qV/l - V 2 \/l - q 2 r? 2 


(27) 


cos' 


V) 


(26) 


11 


normal stress a zz becomes 

W P " + eq3) (28) 

where F is the total force. This result is identical to that obtained by Rae in 
reference 11 . 

Dirac delta at a point . - From equation (27) it is possible to obtain the incident nor- 
mal stress at the wave front due to a Dirac delta loading at a point which is written as 

4Fq3 c 

( a zz)p - -o-H 5 + 2q2 + 16t > 3 ) (29) 

TTZ 0 

where F is the impulse of the load. 

Case of a distributed load where a ^ 0 .- For the case of a distributed load, the 
solution near the wave front z * Cjt will always satisfy the relation 

c^v(^ (30) 

This makes the term under the radical in equation (26) negative and makes rj pure 
imaginary. (In order to consider the stress at any point other than at the wave front, 

77 would, in general, have to be considered a complex number, dependent upon z, a, 
and ip.) Therefore, by letting 77 = iy and rj = -iy where 

(31) 


Cjta cos x/y + zyz 2 + a 2 cos 2 i^ - Cj 2 t 2 
z 2 + a 2 cos 2 )// 


equation (24) becomes 



a 2 P d 2 C 71 

irC]% dt 2 ^0 



(32) 


Near the wave front, z « Cjt and y « 1. Using this knowledge and the fact that 
q < 1 allows equation (32) to be expanded, using the binomial expansion, to yield 

(l - 4q 2 y 2 + 4q 4 y 4 ) + 4q 3 y 2 ^l - &L . J&1+ . • -)(l " ^ ^ + * * *) 

(33) 

12 


(a ) C 

^ ZZ ) P J Vcj 2 dt 2 Jo 



The expansion was truncated at powers of y 2 to yield 



a 2 P d 2 p 

TTCj 2 dt 2 ^0 


yy(4q3y 2 - l) 


sin 2 i^ di (/ 


( 34 ) 


It can be shown from equation (34) that the step load for a = 0 given by equation (28) 
can still be obtained, but the Dirac delta load for a = 0 becomes 


ttFq 3 ^ 


(35) 


which is less than that determined in equation (29). Since the series expansion has been 
truncated, the solution presented in equation (35) is an approximation. Consequently, the 
results of the investigation for the Dirac delta loading is limited to the case where 
a > 0. 

Substituting for y and y from equation (31) into equation (34) and making the 
substitution u = cos i f/ results in an integrand with 12 identifiable terms. Equation (34) 
can be written as 


12 


. \ _ a 2 P d 2 V T 

t zz )p~irc 2 dt 2 l n 


n=l 


(36) 


where I n are the integrals of these 12 terms. A representative example of one of these 
integrals is Ij, written as 




- a 2 u 2 ) (l - u 2 ) ^ du 




2 


(37) 


and the remaining 11 integrals are presented in appendix C. 

In order to evaluate these integrals, an expansion was utilized. Since | u | i 1 
over the open interval -1 < u < 1, then 


(1-u 2 ) 1/2 =1-^-^-... (38) 

This expansion diverges at the end points. However, the end points of the interval in 
the 12 integrals may be excluded since their integrands at the end points remain finite. 
For ease of computation, this series expansion was truncated such that the numerator of 
the integrals always contained terms of u to the fourth power. The expressions for the 
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normal stress (cr Z z)p at tiie wave front for the step load and Dirac delta load are 
shown in appendix D. 


RESULTS 

Normal impulsive loading of an elastic plate by a uniform and point load was inves- 
tigated. Results were obtained for the incident normal dilatation stress wave along the 
axis at the wave front. The time-dependent step and Dirac delta loadings were 
considered. 

Figure 3 shows the typical results of the incident normal stress at the wave front 
for a step load applied uniformly over a circular area as a function of the radius of 
loading a. The total load applied in each case was a constant 27.95 N (2 tt lb), and the 
distance traveled was z = 1.27 cm (0.5 in.). The plate was assumed to be aluminum. 

It will be noted that the solution for a 0, as given by equation (Dl), limits to the solu- 
tion for an a = 0, as given by equation (28). The decrease in normal stress for 
increasing a is rapid when a is small. 

Figure 4 shows corresponding results for a Dirac delta load uniformly applied 

16tQ7T \ 

over a circular area. The load was determined from P = 44.718tQ.N-s IF = — lb-sl, 

where tg is chosen as 2.45 x 10~® s from reference 8 (according to the Hertz theory 
of impact). Again as a is increased, the stress level at the wave front decreases. It 
will again be noted that the solution for a ^ 0, as given by equation (D8), is consistent 
With the solution for a = 0, as given by equation (29). 

Figure 5 shows the nondimensional ratio of normal stress at the dilatation wave 
front to applied pressure, as a function of a/z for an applied step load. This was 
obtained from equation (Dl). The curve in figure 5 shows (as expected) that either as a 
gets large (plane wave) or as z gets small (the pressure application point) the ratio of 
normal stress to applied pressure approaches unity. Another observation that can be 
made is that, for a given z, the plane-wave case (i.e., where no attenuation at the wave 
front is experienced) is approached very rapidly for values of a > z. 

Figure 6 shows the progression of the wave front as it propagates into the plate for 
a step input. The total load applied is again a constant 27.95 N (271 lb) but the intensity of 
the pressure P is varied. At z = 0 the normal stress o zz equals the applied pres- 
sure for a particular radius a. As a increases, the rate of decay decreases, and the 
solution rapidly approaches that of a plane wave. 
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Radius of loading, a, in. 

Figure k.- Incident normal stress along the axis at the wave front for 
Dirac delta input. F = 44.7l8tQ N-s 
(1.0 in.); t = 2.45 x l<r 6 




Z = CEl 


S. 



1.00 p- 



Figure 5*- Normal stress as a function of a/z for a constant pressure amplitude P. 
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CONCLUDING REMARKS 

The normal stress due to the incident dilatational wave produced in an infinite 
plate when loaded normal to its surface has been evaluated. The loading was distributed 
uniformly over a circular area. Solutions have been obtained on the axis at the wave 
front for step and Dirac delta time variations. 

Although the solutions presented are limited, they do show the effect of increasing 
the area of loading. In a similar manner, other contributions to the normal stress, such 
as the incident shear and reflected shear and dilatational waves, could be obtained. How- 
ever, the superposition of one wave upon another cannot be determined by this solution 
since the different wave fronts pass the same point at different times. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 20, 1972. 
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APPENDIX A 

DEVELOPMENT OF TRANSFORMED NORMAL STRESS 

By substituting equations (3) and (4) into equations (1) and (2) and taking the Laplace 
transform, two partial differential equations result, one of which (eq. (14)) is 


8 2 A 1 8A _ _p^_£ 8 2 A 

r 8r 


= 0 


8r c 1 or dz* 

where p is the Laplace transform variable. Separation of variables by assuming 


(Al) 


A = R(r) Z(z) 


yields 


= ^ [ A (|) J 0 (r|) sinh az + B(ij) Jg(r£) cosh azjd£ 


(A2) 


where | 2 = a 2 - -2—, a is a constant, and A(£) and B(|) are arbitrary functions. 
C 1 

Similarly, the other partial differential equation (eq. (13)), when solved, results in 
fi = [c(0 Jj(r£) sinh /3z + D(£) Jj(r|) cosh /3zjd£ 


(A3) 


where! /3 2 = £ 2 + -2—, /3 is a constant, and C(£) and D(£) are arbitrary functions. 


c 2 
c 2 


Substitution of the values for A and ft into the Laplace transform of the differential 
equations (1) and (2) yields 

p 2 u r = J |c 2 2 £a(£) sinh az - B(£) cosh azj 

+ C2 2 j3|"^|c(^) cosh /3z + D(|) sinh J^(r£) d£ 


(A4) 
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APPENDIX A - Continued 


and 


p 2 u z = cosh az + B(£) sinh azj 

- c 2 2 |c(£) sinh /3z + D(£) cosh /3zj| J Q (r|) d| 
From equations (7) and (8) the transformed stresses are found to be 
a zz = X ^ [ A (0 sinh az + B(2j) cosh azj Jq(t£) d£ 

+ ^2 ^ ^ 2 c^ 2 ^“^|a(|) sinh az + B(£) cosh azj 

- /3c 2 2 [c( 4) cosh /3z + D(|) sinh 0zjj Jg(r£) I d £ 


(A5) 


(A6) 


and 


^rz ~^2 ^ ^ c 1 2 «|a( 4) cosh az + B(|) sinh azj 

+ c^r 1 ^ 2 + ^^|c(4) sinh /3z + D(£) cosh 0zjj J^) £ d£ (A7) 

By applying the transformed boundary conditions, equations (9), (10), and (lib), to equa- 
tions (A6) and (A7) the following equations are obtained: 


-2c 1 2 oAU) + 2c 2 2 r 1 ^ 2 + ^J D (I) = 0 

-2cj 2 a|A(^) cosh ah + B(|) sinh azj 
+ 2c 2 2 r 1 ^ 2 + ^^[CU) sinh 0h + D(|) cosh /3hj = 0 


(A8) 


(A9) 


X + 2g a 2 c 2 W) - ^ )34 c 2 2 C(?) = PaJ^) 
k P Z ! P' 2 


(A10) 
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APPENDIX A - Concluded 



^ a 2 Cj 2 j £ _1 |a(£) sinh ah + B(£) cosh 



_ 

n2 


jSc 2 2 |c(£) cosh |3h + D(|) sinh |Sh = 0 


(All) 


Equation (A10) is obtained with the help of the Hankel inversion theorem (ref. 16) and an 
expression from reference 14. Equations (A8) to (All) are solved simultaneously for 
A(£), B(|), C(£), and D(|) which are then substituted into the expression for the nor- 
mal stress, equation (A6). This gives the relation for the transformed normal stress as 


zz 



M(z,£) Jj(a^) J 0 (r£)d£ 
G(i> 


(A12) 


where 


M(z,£) = (cosh az + cosh /3 z)|n.|N 2 (i - cosh ah cosh /3h)j 
+ ^Nj 2 cosh az + N 2 2 cosh /3z)(sinh ah sinh /3h) 

+ (Nj sinh az - Ng sinh 0 z)(N 2 sinh ah cosh /3h - Nj cosh ah sinh /3h) (A13) 


G(£) = (n x 2 + N 2 2 )(sinh ah sinh 0h) + 2N^N 2 (l - cosh ah cosh /3h) 


(A14) 


and 


2mc, 2 


Ni = 

1 r>2 


n2 


4 2 w^ 


2 c 2 2 , 


N 2 = ! .2 


2 mc 1 2 \_„.2 




(A15) 


(A 16) 
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APPENDIX B 


TRANSFORMED STRESS, SHOWING WAVE COMPONENTS 


The transformed normal stress is given by equation (15) as 


^zz = Pa 


,0 ° M(z,£) Jj(a^) J 0 (r|) d£ 

0 GU) 


h = tanh T 1 


t 2 = tanh 


*4 

t 2 


cosh ah = 


sinh ah = 


‘-•i* 




M(z,£) 1 


G(5) ^Nj 2 * N 2 2 j x - NjNj.j'l + > : 2 )J ^ 

+ Ng cosh jSzj^xNg - N^l + j sinh 

-N 2 sinhpz]x N2 (l + t 2 2)-^(l + tl 2) 


2 j ^Nj cosh az|^N 2 (l + x 2 ) 
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APPENDIX B - Continued 
Substituting equation (B4) into equation (Bl) gives 

°zz = § [jl N l cos ^ az + ^2 N 2 cosl1 + ^3( N 1 sinl1 

- N 2 sinh *)] J Q (r|) J^a^) d£ 


where 


az 


Nr 


F i'n 1 + n 1 (n 1 x-n 2 ) + n 1 -*n 2 


(B5) 


(B6) 


F 1 , N 1 T , 1 

2 n 2 N 2 (N 2 x - N t ) N 2 - X N X 


(B7) 


Fn = — — — r + 1 


3"t 


2 (N 2 - X N 1 ) xN 2 - Nj 


(B8) 


Let X = e and Y = e - ^, then 

l + x 


t 2 = — 

2 1 + Y 


-mm) 


(B9) 


Substituting equation (B9) into equation (B6), after some algebraic manipulation, results 
in 


(BIO) 


1 

1 i 

1 

Ni 

N 2 
1 - -s.\ 

f!-X\ 


N 1 

V + x/ 


_Y 

Y 


N l/ l - X \ ( l + Y \ 
' Nn \1 + X/ \1 - y) 
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which, after using the binomial theorem for 


and 

1 + X 


•, becomes 

1 - Y 



(Bll) 


(B12) 


(B13) 


-a z 


Equation (B5) can now be rewritten as 

= f C t F l Nl + F 3 N l) e “ Z + ( F 1 N 1 - F S N l) e ' 

+ ( f 2 N 2 - F 3 N 2) eSZ + ( F 2 N 2 + F 3 N 2) e '^ J 0 < r » d « 
The coefficient of e az , can be expressed as 


1”1 + 3 1 



+ 4XY 




Nj N 2 IN 2 




+ • . • 


(B14) 
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or 


F 1 N 1 + F 3 N 1 = _4X 


K( N 1 + N 

2) 

+ 4XY 

" 2NjN 2 


u 

[ N 1 - N 2) 

l 2 _ 

( N 1 - N 2) 

I 2 


(B15) 


Similarly the coefficient of e~ az is found to be 

PN r N 4Nl I 4X%l(Nl + Nz)2 , 8 y2n i n 2( n » I n 2 ) 

1 1 3 1 N, - N 2 ( Nl - N 2 ) 3 (Nj - N 2 )3 

( N i - N z ) 3 


(B16) 


Thus, from equations (B14), (B15), and (B16) it follows that the stress can be written as 

5 


where 


5 _ = aP 
zz 


Q 0 + Q 0 + Yj ( Qn + Qn *) + • • • 


n=l 




(B17) 


(B18) 




r \ 00 

' N l( N l + N 2 I 

^0 

_( N 1 - N 2) 2 _ 


2 Nl N 2 


( N 1 - N 2 )‘ 


^3 = In 


r \ °0 

N i( N i * n 2 ) 2 

^0 

_ ( N 1 - N 2) 3 _ 


J 0 (r5)J 1 (a?)e'“ (2h ' z) d{ 


J 0 (r{)J 1 (a{)e' h(ai+ ^“ z d| 


J 0 (r|) Jj(a{) e- a < 2h+:z ) d{ 


v-i; 


2N t (N 2 2 + 3N t N 2 ) 

( N l- N 2) 3 


(B19) 


(B20) 


(B21) 


J 0 (r£) Jj(a^) e Ma+0)-az (B22) 
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2 Nl N 2 (Ni + N 2 ) 
( N l- N 2) 3 


J 0 (r9 Jj(a|) e' 2, ' Jh ‘“ z 




(B23) 


The values of Q * (n = 0,1,2, 3,4, 5) are obtained from Q n by interchanging 

with N 2 and a with /3. The results given by equations (B17) to (B23) are identical 

to those presented by Thiruvenkatachar in reference 12. 
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APPENDIX C 


TWELVE INTEGRALS USED IN EQUATION (36) 


In equation (36) the solution of the normal stress was shown to be 


1 n=l 

The twelve integrals of this equation are as follows: 

^1 c^t^z 2 - a 2 u 2 )(l - u 2 ) ^ du 
' l (z* + aV) Z 

p c 1 ^zat^u(l - u 2 ) ^ 

2 ’ •'- 1 (z 2 + a 2 u 2 ) 2 (z 2 - Cj 2 t 2 + a 2 u 2 ) 1/2 

pi c^zau(l - u 2 )^ 2 ^ 2 - Cj 2 t 2 + a 2 u 2 J 
' 3 (z* + a 2 u 2 ) 2 


’I 


1 4q 2 Cj4t 2 a2u 2 (a 2 u 2 - z 2 )(l - u 2 )*^ du 


(z 2 + a 2 u 2 ) 


pi 4q 2 z 2 c (a.2 u 2 _ z 2 )^z 2 - + a 2 u 2 j (l - u 2 ) ^ du 

15 = ^ ( Z 2 + aV) 4 


pi 8q 2 c^ 2 t 2 zau(z 2 - a 2 u 2 ) ^z 2 - Cj 2 t 2 + a 2 u 2 ) ^ (l - u 2 ) ^ du 

^ l S^wj 1 (C7) 




1 -4q 2 Cj5t4 a 3 zu 3 ^ _ u 2)*/ 2 <} u 
^ (z 2 + a^u^) - c^2t2 + a2u2^' 
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APPENDIX D 


EXPRESSIONS FOR THE NORMAL STRESS 


The normal stress on the axis at the wave front z = Cjt for a step load is given 

by 



a 2 P d V T 

,TC 2 dt L “ 
1 n=l 


(Dl) 


where, at z = Cjt 




z 2 ) (4a 2 + 3z 2 ) 
z 2 + a 2 


tan-l(a/ Z )r 2 ( 3z 2 _ 16a 2 

az L 



(D2) 



6^*0 ^z 2 

z 2 (z 2 + 4a 2 ) 

1 {lz 2 + 16a 2 ) j 1 ( tan _ *(a/z) 

a 2 

j(z 2 + a2) 3 

12(z2 + a 2 ) 2 8(z 2 + a 2 ) 8az 


(D3) 


dig q 2 z 2 c 2 
""dt" = 


1 p-f- i T 3 z 2 (l6a4 + 8a 2 z 2 - z 4 ) + z 4 (z 2 - 4a 2 ) 

2 * 4 \3(z 2 + a 2) 3L J 

=-[3z 2 (32a 4 - 32a 2 z 2 + 7z 4 ) + z 2 (-7z 4 + 4a 2 z 2 - 96a 4 jl 

12z 2 (z2 + a 2 ) 2L J 

8z 4( z 2 +a 2) |! z2(32a4 - + z2 t z2 + 8a2)(z2 ' ^ 

t _ an-l(a/z) ^ 2 ( 3 2a4 _ z 4) + z 2( z 2 + 8a 2)( z 2 . 4a 2)j\ 


(D4) 


dljj -12q^Cj 2 z 4 
“d ~ = a 2 


2a 2 + z 2 + 2a 2 + 7z 2 + 2a 2 - z 2 + tan~ 1 (a/z) /o a 2 _ z 2j 

3(z 2 + a 2 ) 3 12z2(z 2 + a 2 ) 2 8z 4 (z2 + P] 8az 5 


(D5) 
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dl^g 4q 3 z 2 Cj2 
~dt ^2 


— -£z 4 + 2z ^( z ^ + 2a 2 )J 


3(z2 + a 2 )‘ 

- -|z 2 (7z 2 - 12a 2 ) - 3z 2 (7z 2 + 2a 2 )] 

12z 2 (z 2 + a 2 ) ^ 

— 2yjz 2 ( 4a 2 - z 2 ) + 3 z 2 (z 2 - 2 a 2 )j 

2( 4a 2 . z 2) + 3 z 2( z 2 . 2a 2)j| 


= ^3- _ S . fa _ f9 . fio 

dt dt dt dt dt dt dt 


= 0 


(D6) 


(D7) 


The normal stress at the wave front z = Cjt, for a Dirac delta loading is given by 

( CTzz )p = ZT2^2 Z In (D8) 


'P tic 2 dt 2 , 
1 n=l 


where, at z = c^t, 
d 2 


I . d2l 2 A - d2l 6 . d2l 7 d2l 8 A d2l 10 

dt 2 1 dt 2 dt 2 dt 2 dt 2 dt 2 dt 2 dt 2 


(D9) 


^2t 


dt* 


12q 3 c 1 3 z 

z 2 (z 2 + 4a 2 ) 

7z 2 + 16a 2 ! 1 j tan _4 (a/z) 

a 2 

3(z 2 + a 2 ) 8 

12(z 2 + a 2 ) 2 8 (z 2 + a 2 ) 8az 


(DIO) 


dr 


3q 3 c 1 3 z 3 

’ . 1 

a 4 

3(z 2 + a 2 ) 


o(^a^ + 8a 2 z 2 - z 4 ) + 


12z 2 (z 2 + a 2 ) 2 

-rA- , v (32a 4 - z 4 ) + tan- 1 (a/z) ( 32a 4 . z 4) 

8z 4 (z 2 + a 2 ) 8az^ 


(32a 4 - 32a 2 z 2 + 7z 4 ) 


(Dll) 
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24q 3 c 1 

a 2 



2a 2 + z 2 t 2a 2 + 7z 2 [ 2a 2 - z 2 

3( Z 2 + a2) 3 12z 2 (z 2 + a 2 ) 2 8 Z 4( z 2 + a 2) 


+ tan-Vz) ( 2a 2 . z 2) 

8az 5 


(D12) 


d 2 I 12 24q 3 Cj 3 z 3 

z 2 + 2a 3 7z 2 + 2a 2 ( z 2 - 2a 2 

dt 2 a 2 

3(z 2 + a 2 ) 3 12z 2 (z 2 + a 2 ) 2 8 Z 4( Z 2 + a 2) 


+ t-HCliaM (-2a 2 + z 2 ) 
8az 5 


(D13) 
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